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Abstract We study properties of quantum chains in a gas of polar bosonic molecules 
confined in a stack of N identical one- and two- dimensional optical lattice layers, 
with molecular dipole moments aligned perpendicularly to the layers. Quantum 
Monte Carlo simulations of a single chain (formed by a single molecule on each 
layer) reveal its quantum roughening transition. The case of finite in-layer density 
of molecules is studied within the framework of the J-current model approxima- 
tion, and it is found that A^-independent molecular superfluid phase can undergo 
a quantum phase transition to a rough chain superfluid. A theorem is proven that 
no superfluidity of chains with length shorter than is possible. The scheme for 
detecting chain formation is proposed. 
PACS numbers: 



1 Introduction 

Ultracold polar molecules, interacting via dipolar forces which can be tuned and 
shaped by external static and microwave electric fields ', represent a fertile ground 
for the study of exotic phases of matter. In recent experiments, high phase space 
density of both polar and non-polar molecules in stable ground states have been 
realized-'-. The JILA group has also reported on observations of dipolar colli- 
sions and an evidence of spatial anisotropy of dipolar interaction in a sample of 
KRb molecules^. The Innsbruck group has produced a high density sample of Cs2 

l:Institute for Theoretical Atomic, Molecular and Optical Physics, Harvard-Smithsonian Center 
of Astrophysics, 
Cambridge, MA 02138, USA 
E-mail: bcapogrosso@cfa.harvard.edu 

2:Department of Engineering Science and Physics, CSI, CUNY, 
Staten Island, NY 10314, USA 
E-mail: Anatoly.Kuklov@csi.cuny.edu 



2 



molecules in an optical latticoi^. The procedure used in— can be generalized to 
heteronuclear (polar) molecules^. 

There exist many predictions of exotic phases formed by polar molecules - 
self-assembled dipolar crystals, dipolar Bose-Einstein condensates, supersolids 

Another interesting proposal concerns polar molecules in optical lattice bi- 
layers with no inter-layer tunneling and induced dipole moments perpendicular to 
the planes^. Such system is expected to display paired-supersolidity, that is, super- 
solid of composite objects formed by two molecules, each belonging to different 
layer. Given successful experimental advances, these phases are very likely to be 
within reach in near future. Here we consider a multilayered system and study the 
formation of quantum flexible chains of dipolar molecules. 

Physical properties of extended quantum objects are of great interest to many 
areas of physics. Quantum strings described by first-quantized Nambu action are 
at the base of modern high energy physics^. In condensed matter, significant ex- 
amples of non-point-like quantum objects include dynamical stripes in high-Tc 
superconductors^^, quantum vortices in superfluids and superconductors, disloca- 
tions in quantum crystals, and chains of polar molecules'- '^. Self-assembly of 
classical dipolar particles, especially in the context of colloids, e.g., ferrofluids, 
has been investigated quite extensively (see Ref. '? and references therein). Re- 
cently, quantum chains of ultracold bosonic polar molecules have been addressed. 
Thermal and quantum fluctuations of chains of polar molecules forming solid 
structures in highly anisotropic traps were studied by Monte Carlo — . In Ref. 
— , the possibility of Bose-condensation of stijf non-interacting chains of polar 
bosonic molecules formed in a stack of pancakes with zero inter-layer hopping 
has been proposed. In this simplified approach the intra-chain dynamics as well as 
molecular exchanges between different chains have been ignored. The fermionic 
case was studied in Ref.-^. 

In the present work we focus on superfluid properties of flexible and interact- 
ing chains of polar molecules in optical lattices. The single chain as well as the 
chain "soup", with allowed molecular exchanges between different chains, have 
been studied. Our results are based on Monte Carlo simulations by the Worm al- 
gorithm— and can be summarized as follows: i) A single chain, formed by one 
molecule per layer, undergoes quantum roughening transition driven by the pres- 
ence of the optical lattice potential within the layers; ii) An ensemble of flexi- 
ble chains can exhibit superfluidity characterized by off-diagonal long-range or- 
der (ODLRO) in the A'-body density matrix only ; iii) No superfluidity of chains 
shorter than the number of layers N can exist, provided that all A' layers are equiv- 
alent to each other; iv) There exist a quantum phase transition from the chain 
superfluid (CSF) to independent molecular superfluids (N-SF) phase. The latter 
is induced by molecular exchanges between different chains. For two dimensional 
(2d) layers this transition is strongly first-order, while for one-dimensional (Id) 
layers it is continuous. 



2 The model 

We consider a stack of N parallel to each other optical lattice layers - either 2d 
planes or Id "cigars", populated by bosonic polar molecules, with dipole moments 
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aligned by an external static electric field perpendicularly to the layers (along z- 
axis) . At low enough temperature the system is well described by the Hamilto- 
nian: 

H'^-J "ia'^./« + 2 ^ yia-jyriianjy, (1) 
(//■), a ia-jy 

where a]^ , cija are creation-annihilation operators of a boson at site and on layer 
a = 1,2,...,N (Latin indexes label sites within a layer, Greek indexes label lay- 
ers); / stands for tunneling amplitude between neighboring sites within each layer; 
riia = aj^aia denotes the onsite density operator obeying the hard-core constraint; 
Via;iy describes the dipolar interaction matrix elements between sites (ia) and 
(jy). This interaction is characterized by a typical strength = /bl, where d, 
stands for the induced dipole moment along the z direction, and denotes the 
distance between two nearest layers. The anisotropic dipole-dipole interaction is 
repulsive within layers and mainly attractive between layers. Simulations have 
been conducted in the situation of zero inter-layer tunneling. 



2. 1 J-current version 

While we use model ([Tl for studying a single chain formed by one molecule per 
layer, the many body case, i.e., finite density n of molecules on each layer, has 
been analyzed within the approach which is a multi-component version of the 
J-current model—: 



ba 



ba 



(2) 



where 'VaJba = Jb{a+i) ^ J&a; Jfoa Stands for an integer vector along bond b be- 
tween two neighboring sites of the D = d + I space-time lattice in the a-th layer 
(each layer is characterized by dimensions Ly = L, V = I, ...,d, and L^z+i = j3, 
with P the inverse temperature); Jjya obeys Kirchhoff 's conservation law V/^J^^^ — 
(absence of interlayer tunneling implies no currents perpendicular to the layers); 
l/K\ ^ J denotes bare superfluid stiffness of molecules in each layer; K2 ^Vd/J 
describes attractive coupling between neighboring layers; the chemical potential 
fj. is the same for all layers. 

Model (|2]l represents a dual formulation of the Bose-Hubbard model, with J^^ 
carrying a meaning of an element of particle world-line in discretized time. Such 
approach provides qualitative insights on phases and phase transitions occurring 
in the original quantum Hamiltonian. We are utilizing it as a first step toward full 
scale quantum ab initio simulations. In this work we consider the simplest version 
- the interactions are cut off to nearest-neighboring sites on nearest neighboring 
layers and all intra-layer terms are ignored. For A'2 = and at /i = this model ex- 
hibits a phase transition from an insulating state to a superfluid one, the latter being 
characterized by macroscopically large closed loops formed by the J-currents^. 

In the present context of A^-layered geometry there are independent super- 
fluids ( that is, N-SF phase) for K2 = 0. In the opposite limit , i.e. at large enough 
K2, it becomes energetically favorable for world-lines in neighboring layers to 
lock in together, so that VaJ;,^ ~ 0. It will be shown that such locking indicates 
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formation of chains of length A'. In other words, there is no superfluid phase of 
chains of length 1 < M < A^. 

In what follows we will be discussing the hard-core version of the J-current 
model, that is, the condition \ Jba\ = 0, 1 is imposed. This implies that the filling 
factor per each layer is < « < 1 and there is always A'^-SF phase at = 0, T = 0. 
The situation can change in the presence of significant interaction between neigh- 
boring sites within the layers induced by the dipole-dipole coupling. Such interac- 
tion can induce solid phases competing with superfluidity. We will consider such 
possibility in the future. 



3 Quantum roughening of single chain 

Let us start by discussing the case of a single chain formed by one dipole per layer. 
We have performed quantum Monte Carlo simulations of the model ([ij for a stack 
of both Id cigars and 2d layers with linear size L. Hard wall conditions have been 
imposed at the space edges ±L/2 of each layer, with one of the dipoles pinned in 
the middle of the layer it belongs to. 

Since a two-body bound state in bilayers is to exist for arbitrarily small dipolar 
interaction strength—, the chain is to be formed at any finite Vd/J, provided L^oo 
and j3 — > °o. Practically, simulations are conducted at finite L and j3, and therefore, 
there is a threshold interaction V^// for chain formation which is strongly depen- 
dent on L. It will be discussed below that such threshold is independent on L when 
the density of molecules is finite. 

We characterize chain roughness by the square of the gyration radius R^: 

" 0=1.2,...,^ 

where stands for the spatial position of the dipole in the a-th layer and (...) im- 
plies full quantum-thermal averaging. When the interaction is turned off, the posi- 
tions of particles are completely uncorrelated and R^^ = R^ = I^, x^ /L!^ ki dl? /M, 
where x, are the lattice sites. There exist two chain phases - rough and smooth 
- which are separated by a quantum phase transition. In the limit of strong inter- 
actions the chain is smooth (stiff), so that the relative displacements of particles 
are of the order of the lattice spacing a, i.e., 7?^ = R^^f^ w a^, and they concen- 
trate around the pinned dipole. As Vd/J decreases, the chain enters the rough state 
where 7?^ diverges logarithmically with A' (see below), but it remains significantly 
smaller than 7?g. The phase transition is characterized by a change in the excitation 
spectrum - from gapped in the stiff regime to sound-like modes propagating along 
the chain length (the z-direction) in the rough phase. 



3.1 Elastic string model 

In order to gain a qualitative understanding, we model the chain as a gaussian 
elastic string subjected to the lattice potential. Let us consider the case of one 
dimensional optical lattices along the x-direction, that is, the displacement of the 
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String is along x, while the sound propagates along z. The quantum roughening is 
controlled by the weight exp(— 5), where the action in units of h reads as: 



1 fP f^^ 
S=- dT dz 
h Jo Jo 



p fdx(z,T)Y , p .{Mz^Y (2n 

' K I — ^ I -(7cos( — X 

a 



2\ dx J 2 ' \ dz , 

' (4) 

where T is the imaginary time; x{z,t) denotes the displacement of the dipole 
belonging the layer z at imaginary time t; p = ni/b~ is the linear mass density 
{m = /2Ja^ is the effective mass in tight binding approximation); U is the depth 
of the optical lattice potential; Vs stands for the speed of sound along the string. In 
the strongly interacting regime Vs can be expressed in terms of the dipolar inter- 
action strength: w Vd/{pb^). 

Rescaling x = a^/{2n), V^T — ^ T (and /3 — > jSVs), we arrive at the canonical 
form of the sine-Gordon model— action: 



SsG^ f dT [ dz- 



J-oo SnK 



dx 



- Ucos(j) } (5) 



where K = is the Luttinger liquid parameter and tj = U / tiVs- The roughening 
transition of Berezinskii-Kosterlitz-Thouless (BKT) type occurs at A" = A";, = 2, 
so that for K > Kc the cos-term can be essentially ignored. In terms of the mi- 
croscopic parameters of the original Hamiltonian J,Vd,a,b:,, we estimate Kc ~ 

^ w 1, where the numerical coefficient of order of unity has been ignored. 

Starting from (|5]l one can evaluate the gyration radius, Eq.lO, in the rough 
regime, where the cos-term can be neglected. Then, the action becomes gaussian 
and Rl = j dx j dz{x^{z,x)) /{b^Nfi) = J d(0 J dq{x^{q,a))) - Jd(oJdq/{(0^ + 
q^) ^ dq/q, where x{q, (o) is the Fourier transform of x(z, t) and the summa- 
tion over Matsubara frequencies ft) is replaced by integration in the limit T=0. The 
resulting divergence is to be cut off at momentum ^ ^ ~ 1 /b^N. This leads to 
the log-type divergence R^= R^Ki XnN °^ \nN. 

At high T, while the chain still exists, the string action (|5]) becomes classical: 
Sct=fi J dzpVf{dx/dz)'^ /2. Accordingly, the gyration radius Q can be evaluated 
as a simple thermal average R^=R^i = J dz{x'^(z)) / {Nb^ ) . In Fourier space, i?^; = 



.JZn^TbiN/{mV,')^N. 



2nmV~ ■'10 q- 



3.2 Quantum Monte Carlo results 

We have performed quantum Monte Carlo simulations of the model ([T]l in the 
r = limit, and have found that the numerical results for are consistent with 
the log- type divergence with respect to in both Id and 2d cases. Fig(T] shows 
the occurrence of the transition at {Vd/J)c^ 1.1 ± 0.1 for a system of 2d layers 
of linear size L = 36 (we have used a cutoff to the nearest neighboring sites for 
b^ = a). At low interaction strength, approaches the expected value for non- 
interacting particles /?g = 216 (at L=36). For large Vd/J, Rg is independent on the 
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Fig. 1 (Color online) as a function ofV^/J in the two-dimensional case. The dipole potential 
has been cut off to nearest neighboring sites. The system sizes are L = 36 and jS = 36, b~ = a. 
Roughening occurs at {Vd/J)c = f=i 1.1 ±0.1. For Vd/J — > 0, approaches the independent 
layers limit i?Q 216 and, for Vj /y > Vf, it becomes 1. 

number of layers and approaches a^, as one would expect for the stiff chain. We 
have also studied the roughening behavior as a function of the cutoff, and observed 
that it is weakly dependent on it. In other words, the dipole-dipole interaction can 
essentially be replaced by the nearest-layer interaction for a single chain. 



4 Multi-chain condensate 

We now turn to the discussion of the case of finite density n of dipoles per layer. 
The existence of a regime where a single chain is rough immediately raises the 
possibility of superfluidity of an ensemble of chains - CSF. It is important to note 
that, in contrast to the single chain case where there is no threshold for V^/J to 
induce binding of dipoles into a chain, molecular exchanges between different 
chains at finite n destroy CSF in favor of the N-SF phase for Vd /J smaller than 
some finite value. 



4. 1 Chain superfluidity order parameter 

In general, ODLRO of chains of length M can be defined via the M-body density 
matrix Dm, in turn defined in terms of M bosonic operators Xffai'^a) placed on M 
adjacent layers: 



£)m(xi,...,xm;x) = (4>"'"(xi,...,xm)*(x,x,...,x)) , 



(6) 
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where (xi , . . . , ) = V^i (xi ) 1/^2 (x2 ) • . . . Vm (xm ) • If Dm is both, short ranged with 
respect to relative distances of the first set of M coordinates Xi , . . . , Xm and long 
ranged with respect to the distance between the center of mass Rem = [xi + ... + 
Xa/]/M and x, and all others Dp are short ranged, then <^ (xi , ...,Xm) describes a 
condensate of chains of length M. Accordingly, the N-SF phase is characterized by 
ODLRO in Dx (xa; x'„) (a = 1 , 2, AO, with all Dm (M < A')trivially long-ranged 
as they can be factorized into products of Di (note though, that Dm will not be 
short ranged with respect to relative distances of the first set of M coordinates). 

One may wonder whether it is possible to realize a condensate of chains of 
length M < N. This can occur quite trivially, if layers are not equivalent to each 
other. As an example, consider the situation where some of the layers, say the 
upper half, are separated by a larger b^, implying negligible binding interaction. 
Thus, while retaining its integrity in the lower half, the chains will partially decay 
into N/2 single-particle superfluids on the upper half. We do not address such 
situations and consider all layers to be equivalent to each other, with the same 
average filling n per site. We prove that under these conditions one can only have 
condensation of N long chains. 

Formation of ODLRO gives rise to classical fields corresponding to broken 
U(l) symmetries. The N-SF phase is characterized by [U(l)]'' broken symme- 
tries, that is, by (quasi-) condensed phases (pa of A'^ U(l)-order parameters. 
Condensate of M-long chains is characterized by the condensed phase <Pm ~ 
2 M<N 9a- For identical layers, though, any choice of M adjacent layers 
where the condensed chains will form is equivalent (we don't consider a possi- 
bility of exotic long-range interactions which can induce pairing between non- 
adjacent layers). Therefore, assuming spatial symmetry between the layers, one 
concludes that there is at least N — I constraints on the phases, which imply that 
all A' phases are condensed, i.e., the N-SF. Alternatively, the spatial symmetry be- 
tween the layers could have been broken spontaneously so that the chains, say, 
M = A'/2 are formed on the upper A'/2 layers and on the lower A'^/2 layers, with 
the plane of zeros separating the two M = A'^/2-condensates. This option, how- 
ever, is excluded by Feynman's argument forbidding the ground state of bosons in 
time-reversal symmetric Hamiltonian to have zerosJ^. Thus, the only option left is 
M = N, which implies the condensation of the field <P = V/'i(xi)i//2(x2)... V/vCx/v). 



4.2 Monte Carlo results 

We have performed Monte Carlo simulations of the model (|2| within a multi-worm 

version of the Worm algorithm- applied to hard-core particles, i.e. = 0, 1. 
The calculated quantities are: the superfluid stiffnesses py, V = l,...,d, and the 

gyration radius G{N). py are determined by fluctuations of total J-currents Z^^' = 

Y.bJba ^ given spatial direction v of the space-time lattice^. In the independent 
molecular superfluid phase for V^/ = 0, there are no correlations between currents, 

that is, {la'^Iat^) ^ Sa.a'- As Vd is increased, correlations appear indicating a non- 
viscous inter-layer drag effect. Upon further increasing V^i beyond some critical 
value, the rough chain superfluid phase occurs, and, in the thermodynamic limit, 
the inter-layer correlator becomes independent of the layer index. This implies 
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Fig. 2 (Color online) Gyration radius G vs K2 normalized by its value Go ~ 1270 at Ki = 0. The 
parameters are = 1, n = 0.5, N = 10, P = 130, L = 130. The transition N-SF to CSF occurs at 
K2 0.42. 



{[la — ^Ji f') I (L'^l^) ^ 0' that is, full locking of superflows in each layer into a 
single flow, so that no super-counter-flows in different layers occurs. 

It is convenient to introduce the Fourier transform of the total J-currents with 
respect to the layer index: 

= -7Ff L exp{iq,a)la , (7) 



where we have used a vector notations for the currents and their transforms; and 
= 2nn^/N, n, = 0,1,2, ...,N ~ 1. We define the superfluid stiffnesses as 

P<^'fe) = ^(/<^'fe)/^^H-^.))o,, (8) 

where {■■■)dn implies averaging over Dp^. For all practical purposes Eq. ^ is 
equivalent to the well known definition in terms of windings^". 

Due to full locking of superflows in each layer into a single flow, the rough 
chain superfluid phase is characterized by p'^'(^;) ^ Sq-.o- In the N-SF phase, in- 
stead, full locking of superflows is lost but correlations between currents p^^\qz) 
persist, resulting in a smooth dependence of p^^\qz,) on q,, with a non-zero min- 
imum at q. = n. The quantum phase transition occurs when p^^\n) — > and 
j9'^^(0) — ^ const 7^ 0, with increasing L and j3. In other words, in the thermody- 
namic limit the N-SF phase coherence lengths ^ along the layers diverge for all 
harmonics q-. In contrast, in the CSF phase ^ becomes finite for all but ^ = 
harmonics. 
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The second observable we have used to characterize the system is the square 
of the gyration radius G{N), defined as follows for the multi-chain ensemble: 

GiN) = ^'£{[Xa-Xo,,]^)D,. (9) 
a.a' 

G{N) describes the lateral spread of the chain A^-body wave-function. Such def- 
inition is a natural generalization of R^, Eq.©, for the single chain case. In the 
NSF phase G{N) = Go = (1 - l/N)Rl w dL^ /12 in full analogy to RI Fig. |2] 
demonstrates that, at a given N (here N=10), the behavior of G{N) as a function 
of the interaction strength K2 is similar to the case of the single chain shown in 
Fig. [D (Nonetheless, as we shall discuss below, unlike the single chain case, in 
the CSF phase chains remain rough even upon increasing the dipolar interaction 
strength A'2.) Beyond the N-SF to CSF transition point (A'2)c, the gyration radius 
G{N) demonstrates the quantum rough chain behavior, i.e. a log divergence with 
the number of layers N (see below). In order to numerically observe such features 
the system size has to be increased so that L > y^G{N) ^ \/lnN, and L> ^ for 

Our numerical results are summarized in Fig. [3] We have chosen the bare 
superfluid stiffness l/Ki =5.0, and considered the case of hard core bosons. 
The main plot of Fig. |3] shows the transition line for the Id (triangles) and 2d 
(squares) systems of linear sizes L = 130, j3 = 130, A'^ = 10 — 40, and L = 24, 
j5 = 24, A' = 16, respectively. Upon increasing the strength of the interlayer attrac- 
tion (terms oc K2 in Eq.(|2]l), chain superfluidity emerges, with a transition point 
which depends on the interlayer molecular density n. The transition lines have 
been determined according to the criterion described above, i.e., p^^\n) — > and 
jo'^' (0) — ^ const 7^ as demonstrated in Fig. 21 

We have found the transition to be first-order in 2d and second order in Id. 
Larger errorbars in 2d (the upper curve on the main panel. Fig. O are due to the 
hysteretic nature of the transition (shown in Fig.|5]). Moreover data for < « < 0. 1 
(and 0.9 < n < 1) cannot be obtained reliably due to strong finite size effects - the 
coherence length ^ at ^, 7^ diverges exponentially as A'2 ^ in full analogy to 
the divergence of the wave function radius of a particle in a bound state of a weak 
potential in 2d. In other words, the single chain size ^ G{N) becomes exponen- 
tially large, so that strong finite size effects are observed (lines in these regions are 
guides to eyes indicating that the transition value of K2 must vanish in the limit of 
a single chain—). 

We have found that, in contrast with the single chain case, there is no rough- 
ening transition as K2 increases - the gyration radius is always diverging with A'^. 
In other words, chains are always rough. However, such divergence is not charac- 
terized by simple ^ InN growth at large A^. Instead, the G{N) ^ InN dependence 
crosses over to a sub-logarithmical behavior for A^ larger than some A^*. This is 
shown in the insets of Fig. [3] where G{N) is plotted as a function of the number 
of layers A^, in Id (right) and 2d (left). Fits to numerical data clearly reveal the 
crossover at A^* < 30 — 40: The InA'^ behavior (fitted by solid lines) for N <N^ be- 
comes G{N) ^ ln'^{N) (dashed lines) for N > N^,, with c w 0.1 in Id, and c w 0.2 
in 2d. Correspondingly, the superfluid stiffness exponentially decreases with A^ as 
the system crosses over to the sub-log behavior (see Fig. |4] below), as if chains 
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Fig. 3 (Color online) Phase diagram within model (|2) in the vs n plane in 2d (L = 2A,N = 
16, = 24, squares), and in Id (L = 130, /3 = 130, TV = 10 - 40, triangles). N-SF (CSF) lies 
below (above) the transition line. Insets: gyration radii G{N), Eq.(|9}, (error bars are smaller than 
symbols) and their fits (lines) in Id ( 208) and 2d ( /?g f« 67). 



were entering some sort of non-superfluid phase due to the caging effect, i.e., 
inter-chain repulsion induced by their shape fluctuations^. 

Fig-Hdemonstrates how the transition points have been determined (here data 
refer to Id case and n = 0.5). The main plot shows p^-^\qz = as a function 
of K2, with the transition point K2 = (^2)0 ~ 0.42 ±0.03 marked by the verti- 
cal arrow. The inset shows the chain superfluid stiffness p''^(^; = 0) for sev- 
eral values of A^. In the CSF we identify two regimes: (^"2)0 < K2 < 0.6, and 
A'2 > 1 . In the latter, the chain binding energy Vb is much larger than the single- 
molecule tunneling /, and tunneling of a whole chain proceeds as a process of 
N-th order in the ratio J /Vb << 1. This results in an exponential suppression 
P''H^; = 0) exp(— AT/13.2). Accordingly, the temperature scale below which 
p'^' {lz = 0) 7^ is also exponentially suppressed. In the region {K2)c < K2 < 0.6 
the stiffness p''^ = 0) exhibits peaks and significant fluctuations characterized 
by large error bars (for AT > 10 ) and non-exponential dependence on AT. The exact 
nature of such fluctuations is an open question. 
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Fig. 4 (Color online) Spatial counterflow superfluid stiffness p''' (g- = vs & in the Id-layers 
case (L = 130, jS = 130, N = 10, 20). The vertical arrow indicates N-SF to CSF transition. Inset: 
the chain spatial superfluid stiffness p'''((?2 = 0). 
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Fig. 5 (Color online) Hysteresis in the CSF stiffness p{q. = 0) (left pannel) and the counterflow 
superfluid stiffness p(q- = n) vs chemical potenial /x in the 2d-layers case. Simulations have 
been conducted for the sizes L = 24, j8 = 24, A' = 8, and for K2 = 1.65. 
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4.3 Detecting chains 

Chains formation is characterized by two qualitative features in absorptive imag- 
ing: the first is vanishing of the single-particle interference pattern (observed after 
the release of independent condensates); the second, a radical change of the struc- 
ture of the density fluctuations. Emergence of these features can be observed as 
Vd I J is increased above the critical value for a given density n. 

The first feature is due to the loss of ODLRO in the single particle density 
matrix D\. The second one is induced by locking of the densities in each layer. 
Indeed, if no chains are present, the columnar absorptive image is characterized 
by density fluctuations Aritot ^ \fNn, where Nn stands for the total number of 
molecules. When chains are formed, fluctuations occur mostly in their center of 
mass position, that is, Autot ^ \/n deep in the superfluid chain phase. 

5 Inter-layer non-viscous drag 

It is useful to present an interpretation of the system within the framework of 
Andreev-Bashkin effect of non- viscous drag between two condensed components^ 
for the 2d case. Such drag has been shown to lead to exotic vortices in two- 
component bosonic system s In the present context there are such com- 
ponents each described by the phase (pa (x) depending on the space coordinates x 
and the layer index a = l,2,...N. The corresponding kinetic energy of the super- 
flows can be written as an extension of Eq.(l) from Ref.— : 



where Pa.a' the effective superfluid stiffnesses, (written in a form which takes 
into account the equivalence of layers as well as their isotropy, i.e. pa.a' is a scalar 
and in general depends on the distance between layers (Z?-|a — a'|). In the case of 
non-interacting layers {Vd = 0), the N-SF phase is characterized by pa,a' = Po^a,a' 
with po denoting intra-layer stiffnesses (the same on each layer). When the dipole 
interaction is turned on, the finite inter-layer drag emerges. This results in a non- 
trivial dependence of pa,a' on {b^\a — ol'\) manifesting a strongly interacting N- 
SF phase. 

In the CSF phase, the full locking of the superflows from different layers im- 
plies that Pa.a' becomes infinitely ranged along z- In other words, pa.a' is inde- 
pendent on the layer indices and is equal to the quantity carrying the meaning 
of the CSF superfluid stiffness p(''(0) presented in Fig. |4] Resorting to the dis- 
crete Fourier transform p{qz) of p(z), the CSF phase is characterized by p(^;) = 
p''^ (0)5fy_ (where = 2nn-/N, n^=0,\,2,...,N — \ for periodic boundary con- 
ditions along z), that is the N-SF to SCF transition is characterized by emergence 
of infinitely-ranged non-viscous drag between layers. 

5.1 The mean field approximation 

Within the mean field description, the relevant interaction term responsible for the 
discussed quantum phase transition is ^ {(p*){\j/)^ +c.c., where {(p) stands for 




(10) 
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the chain order parameter and ( iff) describes molecular condensates in each layer. 
One can therefore write down the minimal Landau free energy in terms of iff and 
<P fields: 

F = ai\\i/\^ + a2\^\^ + bi\xi/\'^ + b2\'P\'^-cu{'P*¥^ + cc), (11) 

where fli 2 , ^1,2, C12 are some phenomenological coefficients. For > 2 this func- 
tional has relevant negative term of order higher than 4, implying discontinuous 
transition between the N-SF state, where 1/ ^ (and trivially <P 7^ 0), and the 
CSF where iff = and <P j^O. This is consistent with our observation of strong 
first-order quantum phase transition in 2d (at T = 0), with the typical hysteretic 
behavior demonstrated in Fig. |5] It is interesting to investigate how the first-order 
transition transforms at finite T in 2d when the vortex proliferation is characterized 
by strong drag effect as decsribed above. 

In D=l-i-l (at T = 0), the mean field description is not applicable. Instead, 
the bosonization approach— i^^*^ describing Berezinskii-Kosterlitz-Thouless sce- 
nario is more appropriate. Accordingly, no hysteresis has been observed in our 
simulations in the case of Id "cigars". 



6 Discussion and summary 

Polarized bosonic molecules confined in a stack of identical one- and two- 
dimensional optical lattice layers form a chain superfluid phase - CSF at large 
enough dipolar interaction. Such phase is characterized by ODLRO only in the A'^- 
body density matrix. Moreover, it is proved that no superfluidity of chains shorter 
than can exist, provided all layers are identical to each other. 

Upon decreasing the interaction strength the system undergoes a quantum 
phase transition from CSF to the A'-independent molecular superfluids. In 2d this 
transition is of strongly first-order (for A^ > 2), while in the Id case, it is continu- 
ous and can be interpreted as a multi-layer version of the Berezinskii-Kosterlitz- 
Thouless transition. 

In terms of an effective description, the quantum phase transition in the 2d 
case is characterized by emergence of the non-viscous inter-layer drag effect with 
infinite range correlations between layers. An interesting question is about the na- 
ture of vortices and their unbinding in the N-SF phase close to the CSF at finite 
T. Similarly to the 2-component cas e^^'^^ , such vortices should be composite ob- 
jects. A closely related question concerns the fate of the first-order transition at 
finite T. 

We have also studied single chain behavior by means of quantum Monte Carlo 
simulations, and have observed the quantum roughening transition from a regime 
of stiff chain, at large dipolar interaction strength, to rough chain, at smaller 
strengths. At zero temperature the rough chain behavior is characterized by log- 
arithmical divergence of the gyration radius with the number of layers A^. In the 
stiff (smooth) phase the gyration radius becomes independent of A'^. 

While no roughening transition was found in the CSF, the single-chain transi- 
tion leaves its "footprint" in the CSF phase: the divergence of the gyration radius 
as a function of A^ crosses over from ^ InN to ^ In'' N type, where c w 0. 1 — 0.2 
(in contrast with the value c = for a single chain). In such a regime chains stay 
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rough, which means that CSF is always a superfluid o{ flexible chains). Crystalline 
order of chains, which would imply loss of roughness, though, is expected to ap- 
pear for strong enough dipolar interaction once intra-layer repulsion is considered. 

In future we would like to investigate the interplay between rough chain su- 
perfluidity, independent molecular superfluids and solid phases by describing the 
system using the full quantum hamiltonian. In particular, the possibility of a su- 
persolid phase in such systems will be explored. Another interesting question con- 
cerns how finite tunneling between layers (which has been neglected in the present 
analysis) can affect the system. 

As far as the practical realization of the chain superfluid phase is concerned, 
we note that our analysis valid for bosonic molecules can be applied to the case of 
fermionic molecules (KRb24) as well with some modifications. Specifically, the 
symmetry between layers should be broken by grouping them into pairs, so that 
the distance between layers within each pair is smaller than the distance between 
the pairs. In this case, it should be possible to achieve strong pairing of fermionic 
molecules from each pair of layers into essentially bosonic di-molecules, which 
in their turn can form chains. 
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